






















microcopy resolution test chart 

NATIONAL BUREAU Of STANDAROS-1963-A 






















































ABSTRACT (Continu'd) 


influence of the neighboring bubbles is relatively short-ranged and that a small number of bubbles can 
adequately model the essential physics in a large array of bubbles. For moderately short separation 
distances, it is found that the E X B rise velocity is substantially reduced in comparison with the 
single-bubble case and that the rise velocity is strongly sheared resulting in the deformation of the 
contours. The implications of the new morphological results on the stability and dynamical behavior 
of the bubbles are discussed. The analysis can also be applied to a multi-plasma density enhancement 
(striatlon fingers and plasma clouds) system such as one might encounter in plasma cloud striation 
fingers. - 
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THE MORPHOLOGY OF A MULTI BUBBLE SYSTEM 
IN THE IONOSPHERE 


I. INTRODUCTION 

The behavior of Ionospheric plasma has been extensively studied with 
the aim of understanding various ionospheric phenomena including equatorial 
spread F (ESF) and plasma cloud strlatlons (see the following reviews: 
Ossakow, 1979; Fejer and Kelley, 1980; Ossakow, 1981; Kelley and McClure , 
1981; Ossakow et al. , 1982). In particular, ESF is thought to be Initiated 
by the Rayleigh-Taylor instability, first proposed by Dungey (1956). In 
the context of this Idea, the plasma density depletions resulting from the 
instability acquire upward polarization Induced E_ x B_ drift velocities. A 
large body of literature has since been developed to describe the linear 
and nonlinear properties of Rayleigh-Taylor plasma density irregularities 
under various assumptions [ Haerendel , 1975; paisley et al. , 1972; 
Chaturvedi and Kaw , 1975a,b; Hudson and Kennel , 1975.]. In particular, 
considerable attention has been given to the morphology and motion 
[ Scannapleco and Ossakow , 1976; Ossakow and Chaturvedi, 1978; Ott, 1978; 
Hudson , 1978; Anderson and Haerendel , 1979; Ossakow et al. , 1979; Zalesak 
and Ossakow, 1980; Zalesak et al., 1982] of plasma depletions ("bubbles") 
in the equatorial ionosphere. Moreover, a number of observations [ Woodman 
and LaHoz, 1976; Kelley et al., 1976; Hanson and Sanatanl, 1971; McClure et 
al. , 1977; Szuszczewicz et al. , 1980, 1981] have indicated the presence of 
rising plasma bubbles. 

Another phenomenon of interest Is that of the strlatlons in 
(artificial) plasma "clouds" (density enhancements). This effect has been 
attributed to the E^ x gradient drift instability [ Llnson and Workman . 
1970; Volk and Haerendel, 1971] and appears to be amenable to treatments 
similar to ESF [Scannapleco and Ossakow , 1976; Scannapleco et al.. 1976; 
Ossakow and Chaturvedi, 1978]. The Rayleigh-Taylor instability and the E_x 
B_ gradient drift instability are both interchange modes that may occur in 
the leading edge of plasma bubbles and backside of clouds respectively. In 
both cases, the resulting density depletions and enhancements are thought 
to drift by the polarization induced E_ x B_ drift. The electric field is 
produced by the polarization of the plasma across the earth's magnetic 
field. The essential ingredient Is the small but nonzero ton-neutral 
collision frequency. 

Manuscript approved November 24,1982. 
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As a result of the Initial Instability such as the Raylelgh-Taylor 
Instability, an array of density depletions are formed with the wave '’Actor 
perpendicular to the earth's magnetic field. Moreover, the leading edge of 
a bubble Itself is Raylelgh-Taylor unstable resulting in further 
"bifurcation" [ Ossakow and Chaturvedl, 1978]. Similarly, the x _B 
instability causes the backside of an Initial cloud to striate, forming 
"finger-llke" structures. Thus, the plasma depletions and enhancements 
typically occur In multitudes. Indeed, McClure et al., (1977) and 
Szuazczewlcz et al. , (1980;1981) have given evidence for multiple bubbles. 
In the previous theoretical and numerical simulation works [ Ossakow and 
Chaturvedl, 1978; Anderson and Haerendel , 1979; Zalesak and Ossakow , 1980; 
Overman et al., 1982], the morphology and evolution have been studied using 
one-bubble models or a uniformly distributed array of bubbles (clouds). As 
a result, the intrinsic influence of neighboring bubbles on each other has 
not been well quantified. Because the bubbles rise due to the polarization 
Induced E_ x B_ drift, the electric field configuration In and around the 
bubbles Is of central Importance. 

In the present paper, we study the structure of electrically 
Interacting multi-bubble systems and seek to Identify the nature and 
effects of the mutual Interaction. For this purpose, we start with a 
simple two-bubble configuration in the context of a fluid description and 
with the emphasis on identifying the basic physics. A simple piecev»is« 
constant density profile is used and the electric field inside and outside 
the bubbles Is obtained analytically by the method of images. The method 
Is then applied to multi-bubble configurations. In a sense, the present 
work is a non-trivial extension of the single bubble work of Ossakow and 
Chaturvedl (1978). 

The scope of the present paper Is limited to discussing the morphology 
of multi-bubble systems and the time-dependent evolution Is not explicitly 
considered. However, the Implications of the results will be discussed In 
the context of the behavior of bubbles consistent with the approximations 
used In the analysis. Although we primarily discuss the ESF plasma density 
depletions (bubbles), the technique developed here can be extended 
straightforwardly to the treatment of plasma density enhancements (clouds). 
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In Section II, we develop a two^bubble model based on an electrostatic 
(dielectric) analogy and solve the current conservation equation with 
quaal-neutrallty for the electric field. Section III describes a multi- 
bubble model and section IV contains the discussion. 

II. A TWO-BUBBLE MODEL 

A. Formulation 

In the present paper, we consider the electric field configuration of 
electrically interacting multi-bubble systems Imbedded In a uniform 
background plasma and neutral gas. In order to illustrate the basic 
physics and the theoretical method, we first develop a simple two-bubble 
model. For his purpose, we adopt a sharp-boundary density profile in which 
the plasma density is piecewise constant, being uniform Inside (n^) and 
outside ( 02 ) the bubbles and having a discontinuity at the bubble 
boundaries [ Haerendel, 1973; Ossakow and Chaturvedt , 1978; Overman et al. , 
1982]. The bubbles are modelled by two-dimensional cylinders at the same 
altitude with circular cross-sections and the axis of the cylinders are 
aligned with the earth's magnetic field which is assumed to be uniform 
along the positive z-axis. No neutral wind Is included. Figure 1 shows 
schematically the geometry and the coordinate system. The bubbles are 
located at x * - x Q and x - x Q so that the inter-bubble separation distance 
is 2x Q . For reference purposes, we denote the cylinders by Cl and C2, 
respectively. The radius of each bubble is a. In the equatorial and low 
latitude regions, the x-axis is along the east-west direction and the 
gravitational force is along the negative y-axis. 

The basic equations in our model are the particle conservation, 
momentum conservation and current conservation equations [see, for example, 
Ossakow, 1981]. In addition, quasi-neutrality can generally be assumed on 
the time scales of interest. In the present paper, we study the morphology 
of two- and multi-bubble systems by solving the current conservation 
equation to obtain the instantaneous electric field perpendicular to the 
magnetic field. 
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The current conservation equation with quasi-neutrality Is 


where £ is the plasma current density due to Ion and electron drifts. 
Neglecting the Inertia terras in the momentum equations for cold Ion and 
electron fluids, the current density In the lab frame can be expressed as 
(see, for example, Ossakow et al. , 1979) 


i’entriXi + ^^i + iE)]. 


where - eB/m^c Is the Ion cyclotron frequency, £ Is the unit vector 
along the earth's magnetic field, £ Is the gravitational acceleration, n is 
the plasma density, and Is the lon-neutral collision frequency. In 
(2), we have used the fact that m e /m^ « 1 and have neglected the 
electron £ x £ contribution. In addition, In arrtvlng at (2), we have made 
the approximation v^/Q^ « l. In the F region, v^/Q^ 1« typically of the 
order of 10“ 2 or less. 

The second term In the square brackets In (2) gives the force- 
field (jj, and E) aligned drift currents due to the finite Lon-neutral 
collisions. It is convenient to separate the electric field according to 


The term ( -a^/e ) £ is the component of E cancelling the drift along the 
gravitational field so that the net drift perpendicular to the magnetic 
field Is described by £. Then, the current J can be written as 








Here, lJ and £ * - g£ are assumed to be uniform. Physically, Ej can be 
thought of as the electric field driving the current £ perpendicular to the 
magnetic field in the frame moving with the velocity Vj * £ * £ 

relative to the lab frame and 0 may be identified as the Pedersen 
conductivity due to * 0. 

Equation (1) can now be written in the equivalent form 

l * (nE) - - ~~ (jxi) • 7n. (5) 

in 

Perkins et al. (1973) obtained this expression and noted that (5) describes 
a dielectric immersed in a uniform electric field E where 



In an earlier work, longmire (1970) utilized a similar magnetostatic 
anology to treat the motion of Isolated ion clouds. Ossakow and Chaturvedl 
(1978) used (5) to analytically study a single bubble system. In this 
dielectric analogy, £ is the polarization (self) electric field of the 
bubbles in the uniform field and E^ corresponds to the total electric 
field (E^j + E) satisfying the boundary conditions across the bubble 
boundaries. 


(oEj), ■ continuous 

(7) 

(Ej)| ■ continuous 
and at infinity (x,y ♦ •) 

Ij + V (8) 

Note that we have implicitly chosen a reference frame in which the electric 
field of the distant undisturbed Ionosphere is E^j. The symbols 
I and 1 refer to the directions parallel and perpendicular to the boundary 
surfaces, respectively. In the present paper, we also adopt the dielectric 
analogy and solve the current conservtlon equation (1) subject to the above 
boundary conditions (7) and (8). For this purpose, it is illuminating to 







rewrite equation (1) as 


V-CoEj) - 0 (9) 

In the following section, we describe the method of Image dipoles used to 
solve this "dielectric equation”. As a matter of notation. In the 
remainder of the paper, we use 1L without the subscript J to denote the 
solution of (9). 

B. The Method of Image Dipoles 

The problem of solving Poisson's equation (9) with multiple 
disconnected boundaries is generally difficult. However, In the case 
treated here with circular cross-sections, the dielectric analogy allows us 
to construct an exact solution. Consider first a single dielectric 
cylinder of radius a centered at x - 0 and a line charge density q located 
at x * b (| b| >a). It Is well known ( Smythe, 1968) that the Induced 
electric field outside the cylinder Is that of a line charge q'» 

- q(l-K)/(l+K) located at x - a 2 /b and a line charge -q' located at x»0. 
The quant tty K Is the ratio of the dielectric constant (o^) Inside the 
cylinder to that outside the cylinder (o 2 ) 



The induced electric field inside the cylinder is that of a single line 
charge q" ■ 2q/(l-HC) located at x ■ b. If we replace the line charge q by 
a line dipole moment - P^x which Is equivalent to two equal and opposite 
line charges separated by a vanishingly small distance, then we find that 
the electric field due to polarization of the cylinder is that of a single 
Image dipole IP given by 

£ “" < ll > 

located at x - a 2 /b. Note that no Image dipole Is present on the axis (to 
be contrasted with the line charge case) and that P^ and £ are collnear, 
pointing In the opposite directions. The Induced electric field Inside the 
cylinder Is that due to a dipole moment £* given by 




located at 


(12) 


■rarl 


We now consider two identical dielectric cylinders of radius a 
centered at x ■ - x Q and x ■ x Q (see figure 1) immersed in a uniform 
electric field E <j . Suppose, for the moment, that the two cylinders are 
non-interacting. Then, for the purpose of calculating the polarization 
electric field outside the cylinders, each cylinder may be replaced by a 
dipole moment P^ located at x « -x Q and x • x Q , where 


(13) 


and K is defined by equation (10). The components of the self electric 
field are 


E x “ 2p 0 t f < x+ V y) + f < x ' x 0 * y) ^ 


(14) 


and 


E y “ 2P 0 [h(x-hc o ,y) + h(x-x o ,y)j. 


(15) 


where 


and 


f(x,y) H - 


(x 2 +y 2 )2 


(16) 


h(x,y) 


(17) 


Here the x- and y-components of the electric field due to a dipole P 

at x - x Q is 2P 0 f(x-x 0 ,y) and 2P Q h(x-x 0 ,y), respectively. The electric 

field E* inside the bubbles is 


E* 


(18) 


In the remainder of the paper, asterisks will be used to denote the 
electric fields inside the bubbles and dipole moments producing the fields. 
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We now allow tite two cylinders to Interact with each other. In 
addition to the uniform external field , each cylinder experiences the 
dipole field of the other. As a result, each cylinder Is further polarized 
giving rise to Image dipoles as given by (11) and (12). Iterating the 
method of image dipoles described above, It Is straightforward to show that 
the total field outside the bubbles Is given by 

E x ’ ‘ E o + n*0 2 P J E(x+ V y) + f < x - x n »y>J* < 19 > 


and 


E , -Jo + t>(*-x 1) .y)|. 


where f and h are defined by (16) and (17). Here, for n * 0, 


P = _ (W) ai 

n ^ t+K'' ^2 n-1* 


( 20 ) 


( 21 ) 

( 22 ) 


and 


x = x 


n o 


(23) 


For n * 0, we have x fl » x Q and P Q is given by (10). 

Similarly, the total electric field Inside each cylinder, say, (C2) 
located at x • x Q , Is found to be 


E ;--w E » + ilic Jo 


(24) 


and 


E ;--4-„io 2p » h ‘* + v’'> • <“> 

where P n and ^ are defined above. For the other bubble (Cl), the field is 
obtained by replacing x 0 with -x Q In the functions f and h. 
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It can be seen from (19)-(23) that the convergence properties of these 
series depend on the parameter s defined by 


1+K 


(2x ) 2 


(26) 


Each successive line dipole moment Is reduced by a factor of s and a 
geometrical factor of order unity from the preceding one. Since x Q > a, we 
have Is| < 1/4 for any K. As a result, the series are generally rapidly 
convergent except for very small center-to-center separation distances 
(x Q ~ a). 

C. Applications 


In order to apply the above dielectric results to the Ionospheric 
bubble and cloud problems In accordance with the dielectric analogy 
(equations (5) and (6)), we identify the dielectric constant o with the 
Pedersen conductivity defined In Section II and recall that E Is given by 
equation (6). Note that the Pedersen conductivity is proportional to the 
plasma density so that 



where n^ and ^ are the plasma densities Inside and outside the bubbles. 
With these Identifications, equations (19), (20), (24) and (25) describe 

the electric field of a two-bubble system In the frame moving with V, - 
-1 

-SI j JX£ relative to the earth. Note that our formalism guarantees that 
the field components obtained above exactly satisfy equation (9) and the 
boundary conditions (7) and (8), as can be verified easily. Thus, the 
solution Is unique. It is also worth noting that all the Image line dipole 
moments are aligned with the x-axls and no higher multipoles such as 
quadrupole moment arise. If the bubbles are at different heights, the 
Interaction can still be expressed as a series of Image dipole moments that 
are parallel to the x-axls but the Images are not Induced on the x-axis 
itself. 

As a result of the polartzatton electric field, the plasma bubbles 
(density depletions) E x B> drift. The drift velocity relative to the 
distant undisturbed Ionosphere where the electric field is Eg Is given 
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by V - c (E* - E ) x B/B 2 . In particular, the single-bubble rise 
velocity Is (using equation 18) 

li" “ (t+|) ^7" * (27) 

In 

For bubbles, K < 1 and is upward. This result has been obtained by 
Ossakow and Chaturvedl [1978]. Because £ is assumed to be uniform, a 
single isolated bubble maintains its ctrcular cross-section as it rises 
with the constant velocity V^. The presence of a second bubble, however, 
modifies the rise velocity significantly. In particular, it is no longer 
uniform and the cross-sections do not remain unchanged. We have 
numerically carried out the summations indicated in (19), (20), (24), and 
(25). The results are depicted in figure 2. The solid lines represent the 
self electric field lines or equivalently the Pedersen current lines while 
the dashed lines correspond to the Instantaneous E x drift 
velocity (Vg) for a two-bubble system (2x Q = 2.5a and K - 0 for 100% 
depletion). This figure shows only one quadrant; the actual system is 
spatially uniform in the z-dlrection and is symmetric about the y-z plane 
and the x-z plane. The solid electric field lines are such that the line 
density is proportional to the field strength. This figure clearly shows 
that the electric field inside the bubble is significantly modified from 
the uniform field of an Isolated bubble (see (18)). As a general remark, 
the interaction vanishes as K ♦ 1 and the separation distance increases to 
infinity. 

In this example, the electric field strength at the point A (x ■ x Q -a, 
y - 0) is approximately 1/3 of that at B (x - x Q +a, y ■ 0). The 

prominently nonuniform electric field inside the bubbles has a number of 
important Implications for the dynamic behavior of the two-bubble system. 
The drift lines (dashed) in figure 2 show that different regions of a 

bubble undergo drift in the east-west (x) direction with respect to the 
undisturbed ionosphere. This horizontal drift can be a significant 
fraction of the vertical rise velocity at some points inside the bubble. 
For example, |E*/(E*-E q )| - 0.55 at r - a_, and 0 - x/8 where r is the 

radial distance from x-x Q , a. is Just inside the bubble surface and 9 is 

measured from the point A (9 - 0). Here, x Q - 1.25 and K - 0 have been 
used (figure 2). The reason for the strong divergence in the field lines 
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A drawing of the electric field lines (solid lines) and the x B_ 
drift lines (dashed lines) showing one quadrant of a two-bubble 
system. Tne system Is symmetric about the x-z plane and y-z 
plane. The separation distance Is 2x Q - 2.5a and K - 0 (100% 
depletion). The points A and B are Inside the bubbles at * ■ x Q -a 
and x - x Q +a, respectively. 






Is that the dipole field due to one bubble opposes the internal field of 
the other bubble. 

Figure 3 shows the two-bubble rise velocity (V y ) at a number of points 
Inside the bubble (C2) relative to the uniform rise velocity of a single 
bubble given by (27). The ratio R 2 = V y /V 1 is plotted for the Interior 
points A and B (figure 2) as a function of x Q /a. We note that the rise 
velocity V y is strongly affected by the neighboring bubble for center-to- 
center separation distances (2x Q ) smaller than 5a to 6a. Moreover, 
comparing V y at A and at B, we see that the E x B drift velocity is 
sheared. This, together with the conclusions of the preceding paragraph, 
implies that an initially circular cross-section would not remain ctrcular 
so that the two-bubble system described here does not correspond to a 
steady-state system. This should be contrasted with single-bubble models 
in which steady-state solutions are possible for piecewise constant density 
profiles. Because a neighboring bubble generally introduces nonuniformity 
in the electric field, it seems difficult to construct a steady-state two- 
bubble system unless x Q /a is large. 

In interpreting the results, we note that our results do not include 
the time-dependence. The above field configuration exists if the two- 
bubble system as described is created at some time (say, t - 0). Thus, it 
would be appropriate as a consistent initial condition for the purpose of 
studying the subsequent time-evolution. In practice, the drift lines shown 
in figure 2 are expected to closely approximate the actual evolution for 
some period of time after t * 0 until the distortion changes the field 
topology significantly. However, the general feature of the field with the 
weakest x-component and hence the slowest rise velocity in the region 
nearest to the neighboring bubble should remain unchanged in time. 


III. A MULTI-BUBBLE MODEL 

In the preceding section, we have discussed in detail a two-bubble 
model. In considering a multi-bubble system that may be applicable to the 
ESP phenomenon, the basic physics and the theoretical treatment remain 
unchanged. However, as the number of bubbles Increases, so does the number 
of image dipoles. For an N-bubble system, the n ttl order expression must 
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x 0 /a 

Pig. 3 The ratio R 2 of the vertical rise velocity V y of a two-bubble 
system to the single-bubble rise velocity (equation (27)) as a 

function of x Q /a. The separation distance Is 2x Q - 2.5a and K - 0 
(100% depletion). The curve A corresponds to the point A (x - x Q - 
a) and the curve B to the point B(x - XQ+a) In figure 2. 
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Include N(N-l) n Image dipoles. Fortunately, the Influence of a bubble 
decreases as the Inverse square of the separation distance. Thus, the 
nearest and the second nearest neighbors are expected to have the dominant 
effects. A consideration of equations (11) and (12) shows that the second 
nearest neighbors have effects of the order of 4 _2 E Q - 0.06 E 0 and that the 
third nearest neighbors have effects of the order of 6 -2 E 0 " 0.03 E Q for a 
given value of a/x Q < 1. This means that the Inter-bubble (and also Inter¬ 
cloud) Interaction Is short-ranged and that the third nearest neighbors and 
beyond have no significant Influence. Thus, only a small number of bubbles 
are necessary to model an N-bubble system with N » l. Note that the 
influence of the bubbles still vanishes at infinity, differing from systems 
satisfying periodic boundary conditions. 

The theoretical treatment described is exactly applicable to any N. 
However, In the remainder of this section, we Include up to a total of 5 
mutually Interacting bubbles. The mathematical manipulations Involved are 
analogous to those of the preceding section, resulting In series 
expressions similar to (19), (20), (24), and (25). As before, only dipole 
moments, not higher multipoles, are Induced. Since no new Insight Is to be 
gained by examining the actual expressions, we give below only the 
results. In figure 4, we sho* the electric field and drift configurations 
of a system with three plasma depletions. The three depletions are again 
modelled by cylinders of radius a, located at x - - 2x Q , 0, and 2x Q . Only 
one quadrant Is shown. The neighboring bubbles are separated by a distance 
2x Q ■ 2.5a, and are 100% depleted (K - 0) as before. The "external” 
electric field is E^ given by (6), and the density profile Is piecewise 
constant. The solid lines represent the polarization electric 
field without . The quantity c (E*- E ) x B /B 2 Is then the 
instantaneous drift velocity relative to the distant undisturbed ionosphere 
and Is represented by the dashed lines. Although the three-bubble system 
is different from a two-bubble system In that the former has a central 
bubble about which the system is symmetric, the general features of the 
field and drift configurations are similar as can be seen by comparing 
figures 2 and 4. That is, the field lines and drift lines inside and 
around the end bubbles of the three-bubble system (figure 4) are similar to 
those of the two-bubble system (figure 2) because the dominant influence 
arises from the nearest neighbors. The central bubble In the three-bubble 
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system is affected by two neighboring bubbles with comparable effects. The 
field is markedly reduced in the region nearest a neighboring bubble (point 
B in figure 4). 


Figure 4 also shows the distortion in the electric field which renders 
the system non-steady-state, as in the two-bubble case. The drift velocity 
(dashed lines) has an east-west (horizontal) component that may be a 
significant fraction of the vertical velocity. For example, 

|E*/(E*-E q )| ■ 0.48 at 9 * n/8 from the point B just Inside the boundary of 
the central bubble and |E*/(E*-E q )| ■ 0.55 at 0 * tt/ 8 from the point x ■ 
x Q -a and y ■ 0 just inside the boundary of the side bubble. Thus, in the 
neighborhood of these points, the bubble elements should have significant 
horizontal drifts. 

In addition to the distortion of bubble contours resulting from the 
non-uniform electric field, the figure also shows that the electric field 
is substantially reduced from that of an isolated single bubble. This fact 
is Illustrated in figure 5 which gives = ^ V y^ V l^ for the P olnts A * B and 
C corresponding to x * 0, x * a and x ■ 2x Q + a, all just InsLde the bubble 
surfaces. Here, Vy is the vertical drift velocity of the three-bubble 
system relative to the distant ionosphere. Comparing figure 5 with figure 
3, we note that the qualitative behaviour of the rise velocity is similar 
in both systems, exhibiting significant reduction from that of a single 
bubble system. However, in the three bubble configuration, the electric 
field inside the central bubble is substantially weaker than the two 
neighboring bubbles. Thus, the central bubble has the slowest rise 
velocity. The line D in figure 5 gives the relative vertical rise velocity 
of a five-bubble system calculated at x * a_. Curve D, to be compared with 
curve B, shows that the influence of the additional bubble on the field 
inside the central bubble is small. The reduction in rise velocities is 
increased as the number of bubbles is Increased. However, the Influence of 
the bubbles beyond the third nearest neighbor is small. In figure 6, we 
show the relative vertical velocity R * v y/ v i as a function of N, the 
number of bubbles, for several values of x Q . For a given separation 
distance, the rise velocity decreases with increasing N and levels off for 
N > 3. Thus, a three-bubble system describes well the basic morphology of 
an N-bubble system. 
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Fig. 5 The ratio Rj of the vertical velocity V y of a three-bubble system 
to the single-bubble rise velocity V^, plotted versus x 0 /a. The 
curves A, B, and C correspond to the points A, B, and C In figure 
4. The curve D (dashed) corresponds to Rj ■ v y/ v i at x » a In a 
five-bubble system, to be compared with the curve B. 
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3 
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Fig. 6 


R " Vy/V^ versus N, the number of bubbles, evaluated Inside 
■x 0 ~a foe the two-bubble case and at x * a for all others. 
0. The separation distances (2x 0 ) are (a) 2x Q - 2.5a, (b) 
3.2a, (c) 2x„ - 4.0a, and (d) 2x„ - 10a. 










Figures 3 and 5 show that the vertical drift velocity of a bubble Is a 
sensitive function of the separation distance (2x 0 ) except for relativly 
large values (2x Q /a >6). In the context of the Raylelgh-Taylor 
instability, this implies that the electric field configuration and the 
drift velocities may depend sensitively on wavelengths (2x Q ). 

As a general remark, we point out that the preceding results derived 
for plasma density depletions for which K < 1 are also applicable to plasma 
density enhancements (clouds) for which K > 1. Calculations for the cloud 
case show that the electric field configurations are qualitatively similar 
to that of the bubble case (figures 2 and 4). In particular, the electric 
field inside the clouds experiences the greatest reduction in the regions 
facing the neighboring clouds. However, the boundary condition 

K(E, ), ■ (E implies that (E, ). Is smaller than (E ), by a 
in i out jl in l out l 

factor of K* 1 for clouds. Thu9, for a given separation distance 2x Q> the 
relative distortion in the electric field lines Inside a cloud is less 
pronounced than in a bubble. 

Finally, in figure 7, we have plotted the vertical drift velocity 
versus K 5 n^/^ for a two-bubble (K < 1) and two-cloud (K > 1) system. 
The velocity is calculated at x - + (x Q - a) (point A in figure 2) and is 
upward for bubbles and downward for clouds. The velocity is normalized to 
Vy (equation (27)) for each value of K. Note that the point with K » 1 
does not exist for each line. For K * 1, the ionosphere is not disturbed 
and there is no bubble or cloud drifting vertically. This is born out by 
the fact that the drift velocity vanishes for a single-bubble (cloud) and 
any multi-bubble (cloud) system. Clearly, the velocities vanish 
differently for different separation distances. Mathematically, the ratio 
of the vertical drift velocity V y - c(E x * - E o )/B 0 to the single-bubble 
(cloud) drift velocity has the limit 

v^“ ■ l 1 - rar* 2f(x + v y) l 

as K * 1. The circled points in figure 7 correspond to the absolute value 
of the quantity which has no physical meaning. 
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I versus. K for two-bubble (K < 1) and two-cloud (K > 
The separation distances (2x Q ) are (a) 2x Q - 2.5a, 
and (c) 2x„ - 10a. 


D 

(b) 








IV. SUMMARY AND DISCUSSION 


In the preceding sections, we have solved the current conservation 

equation (1) with quasl-neutrallty for two-bubble and three-bubble systems 
using a dielectric analogy. These two configurations Include the dominant 
near-neighbor Interaction and can model the essential morphology of the 
multi-bubble systems described. This Is demonstrated by actually 
calculating the field Including up to five bubbles. The solutions are 
exact, satisfying the specified boundary conditions on all the multiple 

disconnected boundary surfaces and at Infinity. 

Equations (19), (20), (24) and (25) give the solution of the two- 

bubble system as a superposition of Image line dipole moments. Similar 
expressions are obtained for three- and N-bubble (cloud) systems (N > 3). 
In all cases, the Interaction Is dominated by the nearest neighbors and is 
sensitive to the separation distance (2x Q ) between bubbles. An important 
result Is that the electric field inside the bubbles is generally highly 
non-uniform so that the multi-bubble and multi-cloud systems are not 
steady-state configurations even with piecewise constant density profiles 

(figures 2 and 4). Moreover, the electric field inside the bubbles is 
significantly weaker than that In a one-bubble system so that the vertical 
drift velocity relative to the undisturbed ionosphere Is slower than the 
one-bubble case (equation (27) and figures 3 and 5). For moderately small 
center-to-center separation distances, 2x Q < 6a, the reduction In the rise 
velocity is substantial. This implies that the electric field 

configuration and the 12 x ji drift velocity may depend sensitively on the 
wavelengths of instabilities causing the initial density fluctuations 
(e.g., the Rayleigh-Taylor instability). In particular, the bubbles with 
smaller wavelength-to-radius ratio would rise more slowly and be distorted 
more strongly. This conclusion may be particularly applicable to the 
Initial linear or early-time non-llnear stages of the evolution. It has 
also been shown that In some regions inside bubbles, the horizontal drift 
velocity may be comparable to the vertical drift velocity with 

|E*/(e*-E q )| as large as 1/2 (Sections II and III). 

A corollary that follows from the non-uniform electric field in the 
bubbles Is that the polarization Induced E x £ drift velocity Is sheared 
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and that a component of the electric field parallel to the local density 
gradient Is developed In the leading edges of the babbles. tor the simple 
piecewise constant density profiles used for our analysis, the density 
gradient Is not well defined. In a more realistic density profile, 
however, some previous model calculations indicate that the electric field 
component parallel to the density gradient has stabilizing Influences on 
the Raylelgh-Taylor instability [ Guzdar et al. , 1982] and the £ x IJ drift 
instability [ Perkins and Poles , 1975; Huba et al. , 1982). lore 
specifically, these calculations show that velocity shear preferentially 
stabilizes the short wavelength modes. In light of these results, we 
suggest that the bifurcation behavior of bubbles and clouds may be 
inhibited by the presence of nearby bubbles (clouds). 

The emphasis of our analysis has primarily been on plasma density 
depletions (bubbles). However, plasma density enhancements (clouds and 
strlation fingers) can also be treated In a similar fashion [ Scannapteco 
and Ossakow , 1976; Scannapteco et al. , 1976; Ossakow and Chaturvedl , 
1978]. If a cloud exists in the equatorial or low latitude region , the 
preceding results are all applicable with the replacement of K > 1 (< < 1 
for bubbles) where K - nj/^* In particular, as tK backside an initial 
cloud begins to bifurcate, the small x /a and "mall N results may be 
applicable. One point to note Is that the electric field Inside a two- and 
three-cloud system with a piecewise constant density profile Is similar to 
that shown In figure 2 and 4 with weaker field strength In the regions 
facing the neighboring clouds. However, the boundary condition K ( E i n )^ ■ 
(E out )i implies that (E tn ) L Is smaller than (E Qut ) L by a factor of K _l for 
clouds. As a result, the distortion In the electric field lines Inside a 
cloud Is less pronounced than In a bubble. 
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OICY ATTN J.P. GILMORE 

COMSAT LABORATORIES 
LINTHICUM ROAD 
CLARKSBURG, HD 20734 
OICY ATTN G. HYDE 

CORNELL UNIVERSITY 

DEPARTMENT OF ELECTRICAL ENGINEERING 
ITHACA, NY 14850 
OICY ATTN D.T. FARLEY, JR. 

ELECTROSPACE SYSTEMS, INC. 

BOX 1359 

RICHARDSON, TX 75080 
OICY ATTN H. LOGSTON 
OICY ATTN SECURITY (PAUL PHILLIPS) 

ESL, INC. 

495 JAVA DRIVE 
SUNNYVALE, CA 94086 
OICY ATTN J. ROBERTS 
OICY ATTN JAMES MARSHALL 

GENERAL ELECTRIC COMPANY 
SPACE DIVISION 
VALLEY FORGE SPACE CENTER 
GODDARD BLVD KING OF PRUSSIA 
P.O. BOX 8555 
PHILADELPHIA, PA 19101 

OICY ATTN M.H. BORTNER SPACE SCI LAR 

CENERAL ELECTRIC COMPANY 
P.O. BOX 1122 
SYRACUSF, NY 13201 
OICY ATTN F. REIBERT 

GENERAL ELECTRIC TECH SERVICES CO., INC 
HMES 

COURT STREET 
SYRACUSE, NY 13201 
OICY ATTN G. HILLMAN 

GENERAL RESEARCH CORPORATION 
SANTA BARBARA DIVISION 
P.O. BOX 6770 
SANTA BARBARA, CA 93111 
OICY ATTN JOHN ISE, JR. 

OICY ATTN JOEL GARBAR1N0 
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GEOPHYSICAL INSTITUTE 
UNIVERSITY OF ALASKA 
FAIRBANKS, AK 99701 

(ALL CLASS ATTN: SECURITY OFFICER) 
01CY ATTN T.N. DAVIS (UNCLASS ONLY) 
01CY A-TN TECHNICAL LIBRARY 
01CY ATTN NEAL BROWN (UNCLASS ONLY) 

GTE SYLVANIA, INC. 

ELECTRONICS SYSTEMS GRP-EASTERN DIV 
77 A STREET 
NEEDHAM, MA 02194 

01CY ATTN MARSHALL CROSS 

HSS, INC. 

2 ALFRED CIRCLE 
BEDFORD, MA 01730 

01CY ATTN DONALD HANSEN 

ILLINOIS, UNIVERSITY OF 
107 COBLE HALL 
150 DAVENPORT HOUSE 
CHAMPAIGN, IL 61820 

(ALL CORRES ATTN DAN MCCLELLAND) 

01CY ATTN K. YEH 

INSTITUTE FOR DEFENSE ANALYSES 
400 ARMY-NAVY DRIVE 
ARLINGTON, VA 22202 
01CY ATTN J.M. AEIN 
01CY ATTN ERNEST BAUER 
01CY ATTN HANS WOLFARD 
01CY ATTN JOEL BENGSTON 

1NTL TEL & TELEGRAPH CORPORATION 
500 WASHINGTON AVENUE 
NUTLEY, NJ 07110 

01CY ATTN TECHNICAL LIBRARY 

JAYCOR 

11011 TORREYANA ROAD 
P.O. BOX 85154 
SAN DIEGO, CA 92138 

01CY ATTN J.L. SPERLING 

JOHNS HOPKINS UNIVERSITY 
APPLIED PHYSICS LABORATORY 
JOHNS HOPKINS ROAD 
LAURAL, MD 20810 

01CY ATTN DOCUMENT LIBRARIAN 
OICY ATTN THOMAS POTEMRA 
01CY ATTN JOHN DASSOULAS 

KAMAN SCIENCES CORP 
P.O. BOX 7463 

COLORADO SPRINGS, CO 80933 
OICY ATTN T. MEAGHER 


KAMAN TEMPO-CENTER FOR ADVANCED STUD: 
816 STATE STREET (P.O DRAWER QQ) 
SANTA BARBARA, CA 93102 
OICY ATTN DASIAC 
OICY ATTN TIM STEPHANS 
OICY ATTN WARREN S. KNAPP 
OICY ATTN WILLIAM MCNAMARA 
OICY ATTN B. GAMBILL 


LINKABIT CORP 
10453 ROSELLE 
SAN DIEGO, CA 92121 
OICY ATTN IRWIN JACOBS 

LOCKHEED MISSILES & SPACE CO., INC 
P.O. BOX 504 
SUNNYVALE, CA 94088 
OICY ATTN DEPT 60-12 
OICY ATTN D.R. CHURCHILL 

LOCKHEED MISSILES & SPACE CO., INC. 
3251 HANOVER STREET 
PALO ALTO, CA 94304 

OICY ATTN MARTIN WALT DEPT 52-12 
OICY ATTN W.L. 1MH0F DEPT 52-12 
OICY ATTN RICHARD G. JOHNSON DEPT 
OICY ATTN J.B. CLADIS DEPT 52-12 

LOCKHEED MISSILE 6 SPACE CO., INC. 
HUNTSVILLE RESEARCH 6 ENGR. CTR- 
4800 BRADFORD DRTVF. 

HUNTSVILLE, AL 35807 
ATTN DALE H. DIVIS 

MARTIN MARIETTA CORP 
ORLANDO DIVISION 
P.O. BOX 5837 
ORLANDO, FL 32805 
OICY ATTN R. HEFFNER 

M.I.T. LINCOLN LABORATORY 
P.O. BOX 73 
LEXINGTON, MA 02173 

OICY ATTN DAVID M. TOWLE 
OICY ATTN P. WALDRON 
OICY ATTN L. I.OUGHLIN 
OICY ATTN D. CLARK 

MCDONNEL DOUCLAS CORPORATION 
5301 BOLSA AVENUE 
HUNTINGTON BEACH, CA 92647 
OICY ATTN N. HARRIS 
OICY ATTN J. MOULE 
OICY ATTN GEORCE MROZ 
OICY ATTN W. OLSON 
OICY ATTN R.W. HALPRIN 
OICY ATTN TECHNICAL LtBRARY SKRVIC? 
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MISSION RESEARCH CORPORATION 
735 STATE STREET 
SANTA BARBARA, CA 93101 
OICY ATTN P. FISCHER 
01CY ATTN W.F. CREVIER 
01CY ATTN STEVEN L. GUTSCHE 
01CY ATTN D. SAPPENFIELD 
01CY ATTN R. BOGUSCH 
01CY ATTN R. HENDRICK 
01CY ATTN RALPH KILB 
01CY ATTN DAVE SOWLE 
01CY ATTN F. FAJEN 
01CY ATTN M. SCHEIBE 
01CY ATTN CONRAD L. LONGMIRE 
01CY ATTN WARREN A. SCHLUETER 

MITRE CORPORATION, THE 
P.O. BOX 208 
BEDFORD, MA 01730 

01CY ATTN JOHN MORGANSTSRN 
01CY ATTN G. HARDING 
01CY ATTN C.E. CALLAHAN 

MITRE CORP 

WESTGATE RESEARCH PARK 
1820 DOLLY MADISON BLVD 
MCLEAN, VA 22101 
OICY ATTN W. HALL 
01CY ATTN W. FOSTER 

PACIFIC-SIERRA RESEARCH CORP. 

1456 CLOVERFIELD BLVD. 

SANTA MONICA, CA 90404 
OICY ATTN E.C. FIELD, JR. 

PENNSYLVANIA STATE UNIVERSITY 
IONOSPHERE RESEARCH LAB 
318 ELECTRICAL ENGINEERING EAST 
UNIVERSITY PARK, PA 16802 
(NO CLASS TO THIS ADDRESS) 

OICY ATTN IONOSPHERIC RESEARCH LAB 

PHOTOMETRICS, INC. 

442 MARRETT ROAD 
LEXINGTON, MA 02173 

OICY ATTN IRVING L. KOFSKY 

PHYSICAL DYNAMICS, INC. 

P.O. BOX 3027 
BELLEVUE, WA 98009 

OICY ATTN E.J. FREMOUW 


PHYSICAL DYNAMICS, INC. 

P.O. BOX 10367 
OAKLAND, CA 94610 
ATTN A. THOMSON 

R & D ASSOCIATES 
P.O. BOX 9695 
MARINA DEL REY, CA 90291 
OICY ATTN FORREST GILMORE 
OICY ATTN BRYAN GABBARD 
OICY ATTN WILLIAM B. WRIGHT, JR. 
OICY ATTN ROBERT F. LELEVIER 
OICY ATTN WILLIAM J. KARZAS 
OICY ATTN H. ORY 
OICY ATTN C. MACDONALD 
OICY ATTN R. TURCO 

RAND CORPORATION, THE 
1700 MAIN STREET 
SANTA MONICA, CA 90406 
OICY ATTN CULLEN CRAIN 
OICY ATTN ED BEDROZIAN 

RAYTHEON CO. 

528 BOSTON POST ROAD 
SUDBURY, MA 01776 

OICY ATTN BARBARA ADAMS 

RIVERSIDE RESEARCH INSTITUTE 
80 WEST END AVENUE 
NEW YORK, NY 10023 

OICY ATTN VINCE TRAPANI 

SCIENCE APPLICATIONS, INC. 

P.O. BOX 2351 
LA JOLLA, CA 92033 

OICY ATTN LEWIS M. LINSON 
OICY ATTN DANIEL A. HAKL1L 
OICY ATTN E. FRIEMAN 
OICY ATTN E.A. STRAKtk 
OICY ATTN CURTIS A. SMITH 
OICY ATTN JACK MCDOUGALL 

SCIENCE APPLICATIONS, INC 
1710 COODRIDGE DR. 

MCLEAN, VA 22102 
ATTN: J. COCKAYNE 
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SRI INTERNATIONAL 
333 RAVENSWOOD AVENUE 
MENLO PARK, CA 94025 

OICY ATTN DONALD NEILSON 
01CY ATTN ALAN BURNS 
OICY ATTN G. SMITH 
OICY ATTN L.L. COBB 
OICY ATTN DAVID A. JOHNSON 
OICY ATTN WALTER G. CHESNUT 
OICY ATTN CHARLES L. RINO 
OICY ATTN WALTER JAYE 
OICY ATTN M. BARON 
OICY ATTN RAY L. LEADABRAND 
OICY ATTN G. CARPENTER 
OICY ATTN G. PRICE 
OICY ATTN J. PETERSON 
OICY ATTN R. HAKE, JR. 

OICY ATTN V. CONZALES 
OICY ATTN D. MCDANIEL 

STEWART RADIANCE LABORATORY 
UTAH STATE UNIVERSITY 
1 DE ANGELO DRIVE 
BEDFORD, MA 01730 
OICY ATTN J. (JLWICK 

TECHNOLOGY INTERNATIONAL CORP 
75 WICGINS AVENUE 
BEDFORD, MA 01730 

OICY ATTN W.P. BOQUIST 

TRW DEFENSE A SPACE SYS GROUP 
ONE SPACE PARK 
REDONDO BEACH, CA 90278 
OICY ATTN R. K. PLEBUCH 
OICY ATTN S. ALTSCHULER 
OICY ATTN D. DEE 


VISIDYNE 

SOUTH BEDFORD STREET 
BURLINGTON, MASS 01803 
OICY ATTN W. REIDY 
OICY ATTN J. CARPENTER 
OICY ATTN C. HUMPHREY 
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NAVAL RESEARCH LABORATORY 
WASHINGTON, D.C. 20375 
DR. P. MANGE - CODE 4101 
DR. E. SZUSZCZEWICZ - CODE 4187 
DR. J. GOODMAN - CODE 4180 
DR. P. RODRIGUEZ - CODE 4187 

A.F. GEOPHYSICS LABORATORY 
L.G. HANSCOM FIELD 
BEDFORD, MA 01730 
DR. T. ELKINS 
DR. W. SWIDER 
MRS. R. SAGALYN 
DR. J.M. FORBES 
DR. T.J. KENESHEA 
DR. J. AARONS 
DR. H. CARLSON 
DR. J. JASPERSE 

CORNELL UNIVERSITY 
ITHACA, NY 14850 
DR. W.E. SWARTZ 
DR. R. SUDAN 
DR. D. FARLEY 
DR. M. KELLEY 

HARVARD UNIVERSITY 
HARVARD SQUARE 
CAMBRIDGE, MA 02138 
DR. M.B. McELROY 
DR. R. LINOZEN 

INSTITUTE FOR DEFENSE ANALYSIS 
400 ARMY/NAVY DRIVE 
ARLINGTON, VA 22202 
DR. E. BAUER 

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
PLASMA FUSION CENTER 
LIBRARY, NW16-262 
CAMBRIDGE, MA 02139 


NASA 
Code 961 

GODDARD SPACE FLIGHT CENTER 
CREENBELT, MD 20771 
DR. R.F. BENSON 
DR. K. MAEDA 
Dr. S. Curtis 
Or. M. Dubla 

NATIONAL TECHNICAL INFORMATION CENTER 
CAMERON STATION 
ALEXANDRIA, VA 22314 
12CY ATTN TC 

COMMANDER 

NAVAL AIR SYSTEMS COMMAND 
DEPARTMENT OF THE NAVY 
WASHINGTON, D.C. 20360 
DR. T. CZUBA 

COMMANDER 

NAVAL OCEAN SYSTEMS CENTER 
SAN DIECO. CA 92152 

MR. R. ROSE - CODE 5321 

NOAA 

DIRECTOR OF SPACE AND ENVIRONMENTAL 
LABORATORY 
BOULDER, CO 80302 
DR. A. CLENN JEAN 
DR. G.W. .ADAMS 
DR. D.N. ANDERSON 
DR. K. DAVIES 
DR. R. F. DONNELLY 

OFFICE OF NAVAL RESEARCH 
800 NORTH QUINCY STREET 
ARLINGTON, VA 22217 
DR. G. JOINER 

PENNSYLVANIA STATE UNIVERSITY 
UNIVERSITY PARK, PA 16802 
UR. J.S. MISBET 
DR. P.R. ROHRBAUCH 

DR. L.A. CARPENTER 

DR. M. LEF. 

DR. R. 0I7ANY 
DR. P. BENNETT 
DR. F. KLF.VANS 
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PRINCETON UNIVERSITY 
PLASMA PHYSICS LABORATORY 
PRINCETON, NJ 08540 
DR. F. PERKINS 

SCIENCE APPLICATIONS, INC. 

1150 PROSPECT PLAZA 
LA JOLLA, CA 92037 
DR. D.A. HAMLIN 
DR. L. LINSON 
DR. E. FRIEMAN 

STANFORD UNIVERSITY 
STANFORD, CA 94305 
DR. P.M. BANKS 

U.S. ARMY ABERDEEN RESEARCH 
AND DEVELOPMENT CENTER 
BALLISTIC RESEARCH LABORATORY 
ABERDEEN, MD 

DR. J. HEIMERL 

UNIVERSITY OF CALIFORNIA, 
BERKELEY 

BERKELEY, CA 94720 
DR. M. HUDSON 

UNIVERSITY OF CALIFORNIA 
LOS ALAMOS SCIENTIFIC LABORATORY 
J-10, M3-664 
LOS ALAMOS, NM 87545 
M. PONGRATZ 
D. SIMONS 

G. BARASCH 
L. DUNCAN 
P. BERNHARDT 

UNIVERSITY OF CALIFORNIA, 

LOS ANGELES 
405 HILLGARD AVENUE 
LOS ANGELES, CA 90024 
DR. F.V. CORONITI 
DR. C. KENNEL 
DR. A.Y. WONG 

UNIVERSITY OF MARYLAND 
COLLEGE PARK, MD 20740 
DR. K. PAPADOPOULOS 
DR. E. OTT 


UNIVERSITY OF PITTSBURGH 
PITTSBURGH, PA 15213 
DR. N, ZABUSKY 
DR. M. BIONDI 
DR. E. OVERMAN 

UTAH STATE UNIVERSITY 
4TH AND 8TH STREETS 
LOGAN, UTAH 84322 
DR. R. HARRIS 
DR. K. BAKER 
DR. R. SCHUNK 












